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^ We use a transport-theory approach to construct the static propagator of a gauge boson in a 

plasma with a general axially- and reflection-symmetric momentum distribution. Non-zero mag- 



netic screening is found if the distribution is anisotropic, confirming the results of a closed-time- 
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path-integral approach. We find that the electric and magnetic screening effects depend on both 
the orientation of the momentum carried by the boson and the orientation of its polarization. In 



present some fairly general conditions on the dependence of these effects on the anisotropy. 
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I. INTRODUCTION 



One of the important problems in non-equilibrium QED and QCD is to study how such 
systems approach equilibrium. In principle one should solve the Schwinger-Dyson equation 
to do this at a purely quantum level, but it is impossible to do this without significant 
approximations For many practical purposes one can use transport equations to study 
equilibration of QED and QCD plasmas with particle collisions taken into account. For 
details we refer the reader to the extensive reviews in the literature 

as. 

However the 

scattering kernel in these transport equations can suffer from severe infrared divergences 
due to ultra-soft photons or gluons exchanged in t-channel. 

It is well known that the collective modes in the QED plasma screen the long-ranged 
Coulomb force jj, . The screened potential then has a Yukawa form in the medium: 

p -m D r 

V(r) ex — — , (1) 

where mjj is the Debye screening mass. As a result the infrared behaviour of the exchanged 
gauge bosons is improved. This Debye screening effect can be generalized to non-equilibrium 
QCD systems and used in the transport equations describing quark-gluon plasmas Q. 
However for magnetic interactions the infrared problems persist. In the case of QED in 
equilibrium the magnetic screening mass vanishes not just at one-loop level but to all orders 
in perturbation theory It also has been argued that in QCD the magnetic screening 
mass should be of the order of g 2 T, and cannot be computed perturbatively Q]. 

The situation is quite different in non-equilibrium systems, at least if they are anisotropic. 
Recently the closed-time-path-integral formalism was used to derive the magnetic screening 
mass at one-loop level j^. This showed that the magnetic screening mass in QED is non-zero 
if the plasma has an anisotropic single-particle momentum distribution, i.e. it depends on 
the direction of the momentum. This would suggest that there is no natural infrared problem 
in such a system. More recently, however, Romatschke and Strickland have used the 
Hard Thermal Loop (HTL) approach to show how nonzero magnetic screening effects can 
arise in anistopic plasmas. Their work shows that there can be antiscreening, in the sense of 
negative screening masses, as well as screening. This feature of out-of-equilibrium systems 



has also been noted by Arnold et al. 
have been studied by Mrowczynski 
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11]. It reflects the instabilities of these systems which 



(For a review of these ideas, see Ref. jl^|.) 
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The instability of these systems implies that, even though some components of the magnetic 
interactions are screened, a simple perturbative treatment is still not adequate. 

It is well known that the polarization tensor of a gauge boson derived from the HTL 



15j. The 



approach at leading order can also be obtained from a transport equation P, 
equivalence of the two approaches has also been generalized to the case of anisotropic but 
homogeneous systems jl6[. In this paper we use the transport-equation approach to derive 
the same expression for the magnetic screening mass as was found in j^. The advantages 
of this approach are that the physical picture is more transparent and its classical character 
is emphasized. We present the general form of the propagator of a gauge boson in an 
anisotropic medium. This propagator is central to the scattering kernel of the transport 



equations [p. The structure of this propagator agrees with that obtained by Romatschke and 
Strickland 10] who assumed a more restricted functional form for the anisotropic distribution 

n 

of particles in the plasma. A minor technical difference from Ref. 10] is that we work in 
Lorentz gauge rather than the temporal axial gauge, but we have checked that our results 
are equivalent to theirs. 

The propagator is more complicated than in the equilibrium case, since it depends on the 
the orientations of the field momentum and polarization relative to the axes of the anisotropy 
of the system. As a result the screening or antiscreening of a transverse magnetic field can 
itself be very anisotropic. Under some fairly general assumptions about the form of the 
anistropv which are likely to cover most cases of physical interest, including those studied 
in Ref. |l0| . we are able to derive conditions on the appearance of screening or antiscreening 
in particular orientations. 

The paper is organized as follows. In section II we present the derivation of the polariza- 
tion tensors in terms of the non-equilibrium distribution function from transport equations. 
In section III we discuss the screening of static fields and, in particular, the form of the 
static polarization tensor for systems with cylindrical symmetry. We summarize our results 
and their implications in section IV. 
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II. DERIVATION OF THE POLARIZATION TENSOR FROM TRANSPORT 
EQUATIONS 



The polarization tensor at one-loop level for a gauge boson in a plasma can be derived 

n n 

from classical transport equations |15( even if the system is out of equilibrium [16j|. As a 
simple example, consider first the case of a gas of electrons interacting with a weak, space- 
time dependent electromagnetic field. The single-particle distribution function /(x, p, t) 
obeys the Boltzmann equation, 

a/( g t P,t) + v • V./(x, p, t) + F(x, t) ■ V p /(x, p, t) = 0, (2) 

where v = p/E(p) is the velocity of a particle with momentum p and 

F(x,i) =e[E(x,t) + vxB(x,t)], (3) 

in terms of the electric field E and the magnetic field B. 

We assume that on the distance- and time-scales of interest the plasma is close to homo- 
geneous and so we can expand /(x, p) around an unperturbed distribution function /o(p) 
as follows: 

/(x,p,t) = / (p) + /i(x,p,t). (4) 

This expansion is similar to the one which is usually used in the derivation of the transport 
coefficient except that /o(p) does not need to be the equilibrium Fermi-Dirac distribution. 
The transport equation (J2J) to first order in f% and the electromagnetic fields is 

^ + v ■ VJi = -F • V p / . (5) 

Solving this for fx yields 

/i(x,p,*) = -eV p / (p)- f tf[E(i',x-v(t-('))+vxB({',x-v(t-i'))] 

J — oo 

POO 

= -eV p /o(p) • / dr[E(t - r, x - rv) + v x B(t - r, x - rv)]. (6) 
Jo 

The induced current density is given by 

j; n ^,t)=eJ^vM X ,p,t), (7) 

where = p^/E(p) = (1, v). Using the above expression for /i(x, p,t), we find that the 
induced current can be expressed in terms of the fields as 

3 ; nd (x) = -e 2 J-^- v, V p / (p) ■ / o °° dr[E(t - r, x - rv) + v x B(t - r, x - rv)]. (8) 



Taking the Fourier transform of this we get 



d3p „ u Mp).[e(^ V xb(4 



(9) 



(2tt) 3 11 vk + ie 

To get this result we have regulated the r integral by inserting a factor of e~ eT with e — > + 
and we have used the fact that the Fourier transform of f£° dre~ eT f(x — vr) is given by 
if(k)/v ■ k + ie. For isotropic systems, the gradient V p /o(p) is proportional to p. The 
second term in Eq. Q vanishes in such cases and so there is no response of the plasma to 
a magnetic field. 

The Fourier components of the the fields can be expressed in terms of the electromagnetic 
potential A lJ, (k) as 



E(k) 



B(ife) = ik x A(ife), 
-ikA (k) + ik A(k). 



(10) 



This allows us to express the current in the form 



jf d (k) 



—ie 



d 3 P V P / (P) 



(2tt) 3 "w-H 
d 3 p 



ie 



-ikA (k) + ik A(k) + i(v • A(fc))k - z(v • k)A(fc)] 



(2VT 



k y% • - V p / • A(k) 
v ■ k + 16 



We can introduce a polarization tensor, defined in terms of the induced current by 

3 ; nd (k) = u^A^k). 



(ii) 



(12) 



The polarization tensor obtained from the induced current is of course only the matter part, 
arising from the gas of electrons. The full tensor also includes a vacuum part which can 
be calculated using the standard field-theoretic methods. After some algebra we get from 
Eq. (jllj) the matter polarization tensor in the Lorentz gauge, 



n 00 (A;) 

n 0i (£;) 

U ij {k) 



d 3 p 
(2^)3 

d 3 p 
(2^)3 

d 3 p 



k-V p / (p) 



k — k • v + ie 
k-V p / (p) 



(27T) 



k — k • v + ie 
/o(p 



(Oij - ViVj) - ; -ViV 



i" 3 



(13) 



E(p)™ ~ l ~ J ' k -k-v + ie' 
If the particles are massless (or the typical energies are high enough that we may neglect 
their masses), the velocity v may be replaced by the unit vector p and the energy E(p) by 

IpI- 



To generalize this result to QCD, one should replace the space-time derivatives in Eq. (J2J) 
by covariant derivatives. However, so long as gA a ^ is small, we can neglect the higher-order 
terms that this introduces. We can then describe the QCD plasma by equations with the 
same form as Eq. (J2J). The resulting contibutions to the gluon polarisation tensor have 
same forms as for the photon except for a trivial factor which counts the relevant degrees 
of freedom. In QCD case one simply needs to replace /(p) in the expression (jl3j) for the 
polarization tensor by a sum of quark, antiquark and gluon terms Q], 2Nf(f q (p) + /g(p)) + 
2N c f g (p). The factor of 2 before Nf is due to the spin degree of freedom. 



III. SCREENING OF STATIC FIELDS IN AN ANISOTROPIC PLASMA 



The long-distance behaviour of the static potential is governed by the low-momentum 
behaviour of the static photon propagator. In the static limit, k = 0, the polarization 
tensor Eq. (pH?|) depends the direction of k only. Denoting the static tensor by n M ^(k), we 
find that its components take the forms 



n 00 (k) 
rF(k) 
ti ij (k) 



d 3 p 
(2tt) 3 

d 3 p 
(2^)3 

d?p 



iPl 



iPl 




(Oij - PiPj) + |p|-~ — PiPj 



IPl 



k ■ p — it 



(2tt) 3 

in a plasma of massless particles. 

In the familiar case of isotropic matter, the static polarisation tensor takes the form 



(14) 
(15) 
(16) 

(17) 



where, in the rest frame of the matter, L^ u = —6^5^ and T^ u = —g^ v + k^k" /k 2 — L^ u . 
For comparison with the results for anisotropic matter, we note that the corresponding 
propagator can be written 



D m {0,k) 



-1 



|k| 2 + m 2 D ' 



D ij (0,k) 



|k| 2 + m 2 



rpij 



a 



(18) 



where a is the gauge-fixing parameter (a = 1 in Feynman gauge) and the Debye screening 



mass is, at the order to which we are working, m 2 D = —Re 



n 



00 



= Re 



IL 



Similarly, the 
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magnetic screening mass is 



IT 



Re 



n 7 



(19) 



but this vanishes in isotropic matter as discussed above. 

In anisotropic matter, the polarization tensor cannot be reduced to just two screening 
masses since it contains, in general, off-diagonal elements and it depends on the direction of 
the momentum carried by the field and on the direction of polarization vector. To examine 
this in more detail we make some specific assumptions about the form of the momentum 
distribution /o(p)> which we believe should be appropriate in the context of a relativistic 
heavy-ion collision. For matter near the collision axis, we assume that the system is sym- 
metric about this axis. We therefore consider a cylindrically symmetric form for /o(p)- Also, 
if we work in the local rest-frame of the matter, there is no net flow and we can assume that 
the distribution is symmetric under reflections, 

/o(-p) = /o(p). (20) 

In this case, as we shall see, the mixed space-time components of the tensor are purely 
imaginary. This means that we can still define a Debye mass as above, although this mass 
now depends on direction. If we use Eq. (|T9"j) to define an averaged magnetic mass from the 
space components of the static tensor (fT^j) . then we find that it reproduces the result for 
the screening mass given in Ref. [^] using the more complicated closed-time-path-integral 
method, 



2f e 2 f d 3 p [ fr2|/(P),ri (- C\2\ V P /(P) 

m * (k) = 2 J [1 + (P " k) + [1 " (P " k) ] k-n 



(21) 



iPl k-p 

Note that the derivation of (}2*Tj) in :9] does not rely on the particular form of the distribution 
assumed here. 

Without loss of generality, we can choose the z-axis along our axis of symmetry and the 
momentum of the field to lie in the xz-plane. We can then express k in the form 

k = k p -k + k z z, (22) 

with k p , k z > 0. In terms of integrals over the longitudinal and radial components of p 
and the angle <fi between p and k in the xy-plane, the components of the static polarization 
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tensor can be written 

t2 r oo 



n 00 (k) 

n te (k) 

ti 0z {k) 
Tl xx {k) 



roc poo p27T , 



e 

- /-OO 



Ppk p cos + p z k z - it 



(2tt) 3 

e 2 
(2tt) 3 

e 2 
(2tt) 3 



rfp 2 / p p G^ 



(ip 2 / ppdpf, 



dp z / 



2tt 



2tt 



2tt 



fc^cos^fg 
J9 p fcp cos0 + p z k z - it 

Ppkp cos (fi + p z k z 



p p COS ( 



it 



-Pz 



1 



\[pI + P 



(p 2 z +p 2 p sin 2 0)/ o 



Pp+P* 



fc_gfccostf + *J£ 



nw < k) = (So* /»*■/. W„ 



p 2 + p 2 L 



rr z (k) = 



(27T) 



^P z / Pp^Pp 



2tt 



H ^ 

Ppkp cos (f) + p z k z 


— «e 


(P 2 Z +P 2 P cos 2 0)/ o 




P 2 p+P 2 z 




kp^ COS (f) + k z 


<9/b 

dp z 


' PpkpCos(p + p z k z 


— ie 


Pp/o 

Pp+P* 




/c p fgcos0 + A;, 


dfo 
dp z 



p sin 



A g2 /-oo /"oo 



PpPz 



pI + pI 



k p ^ cos <j> + k z ^ 



Ppkp cos + p z k z - ie 



COS( 



(23) 



The expressions for the Oy, xy and components have been omitted because their inte- 
grands are odd functions of and so they integrate to zero. 

In general, with k p ^ and k z ^ 0, we can define 7 = k z /k p and carry out the integration 
over to get 

-.2 roo 



n 



00 



yoo poo 
2 / ^ / ^Pp' 

dp 



Pp 



P^ +Pp 



/o 



Pz 



Pz+Po 



+ 



n 0z = 



2tt 2 

''2tt 2 y 



dpp- 



Pz 



dfo dfo 

' dp z dpp 

dfo df 

dpz dpp 



k z Pz 



^k 2 z p 2 z -k 2 p p 2 p ' 



poo 

dp z j dp p 

J"fPz 



dfo dfo 

Pp—~Pz — 



k z Pz 



dpz Z dp p \ Jkffi-kffi' 
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n 0x 



-7 n o z , 

o2 r oo 



dp z / dp. 



l 2 Pz 



2vr 2 Jo ^ z Jo P Jpj+ pj 



Pp 



dfo df 



dp z 



Pz 



+ 



2tt 2 Jo 



dp z 



7Pz 



dpc 



l 2 Pz 



pI+pI 



dp PJ 
dfo 



dfo 

Ppa P'a 

op z op p 



k z Pz 



n 



yy 



2tt 2 



dp z 



dp p 



l 2 Pz 



p 2 z + pI I °p 



Pp: 



dfo 



Pz 



dfo 
dpp 



dpz / dp. 



l 2 Pz 



2vr 2 Jo '"Jo *\ Ip2 + pj 



dfo df 



Ppa 



fc 2 p 2 - A;2p2 



A; 2 p 2 - /c 2 p 2 



n, 



2tt 2 Jo 

„2 



dpz 



A P z 

<fpj+pi 



Pp 



dfo df 



+ 



2?r 2 Jo 



dp z 



7Pz 



dp. 



Pz 



pI+pI 



dp z 
dfo 



Pz 



Pp 



dp z 



Pz 



dp pi 
dfo 



k z Pz 



k z Pz 



d Pp \ ^hlpl - kffi 



2tt 2 Jo 

„2 



dpz 



oo 

dpp 



Pp 



dfo 

dp z 



Pz 



dfo 
dp P 



2vr 2 



dp z 



■ypz 



dPt 



IPz 



Pp 



dfo_ dfo 

dp z z dp p 



k z Pz 



k zP 2 z - k pP 2 p 



(24) 



'p;+Pp 

Note that if the distribution function satisfies the reflection condition (}2Tj|) then the time- 
time and space-space components in the static polarization tensor are real and the mixed 
space-time components elements are purely imaginary. (This can be seen by making the 
change of variables p z — ► — p z , 4> — ► 4> + rr in (|23|)). 

The elements of this tensor do not have definite signs and so there can be antiscreening in 
some directions as well as screening in others. Also, the appearance of a nonzero space-time 
component means that in general the Debye screening of electric fields cannot be separated 
from the effects on magentic fields. The structure of this tensor agrees with that obtained 
by Romatschke and Strickland in the temporal axial gauge. However the Lorentz gauge 
used here is more convenient for handling static electric fields. 

For our chosen form of distribution, the spatial part of static polarization tensor is real 
and symmetric. Hence the tensor H 1 - 7 has three real, orthogonal eigenvectors. The condition 
that the polarization tensor be transverse reduces to 



k,TT j = 0, 



(25) 



in the static limit. This shows that k is an eigenvector with zero eigenvalue. By inspection 
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of Eq. (|23|) we see that the ^-direction, perpendicular to the plane containing k and the 
collision axis, is also an eigenvector. Hence three unit eigenvectors for the static tensor are 

ei = k, e 2 = y, e 3 = k x y. (26) 

We denote the corresponding eigenvalues of the spatial part of the tensor by Hi = 0, H 2 and 
n 3 . In the coordinate system defined by these axes the spatial part of the tensor is diagonal. 
Also we find that n 01 = n 02 = 0. The corresponding static propagator for the gauge boson 
has the non-zero elements 

/ 



D 00 (0,k), £> 03 (0,k) -1 



£> 30 (o,k), £> 33 (o,k) j (|k| 2 + n + )(|k| 2 + n. 



|k| 2 + n 3 , 
n 03 , 




DL1 ^ = ^ D22 ^ = WTiW- (27) 

Here we have introduced the two eigenvalues ll-i- of the full tensor in the time-3 subspace. 
These are given by 

1 ^ ^ , . 1 



n± = --(n 00 - n 3 ) ± -V(n 00 + n 3 ) 2 + 4|n 03 | 2 . (28) 

These results show that the screening or antiscreening of static electromagnetic fields in 
an anisotropic plasma can depend on both the orientation of k, the momentum carried by 
the field, and the orientation of the field itself. It is worth considering two special cases. 
First, when the momentum carried by the field is perpendicular to the collision axis (k z = 
and hence e 3 = z) we have 

n! = n 2 = o, n 3 = n-, n 03 = o. (29) 

In this case fields in the ^/-direction are unscreened. In contrast when the momentum carried 
by the field lies along the axis (k p = and hence e 3 = x) we have 

n 1 = o, n 2 = n 3 = n^, n 03 = o. (30) 

In this case the screening masses for the transverse components of the field are degenerate. 
For both of these special directions of the momentum, there is no mixed time-3 component 
and so — n o and H 3 are both eigenvalues. 

In general there can be antiscreening as well as screening of electromagnetic fields in an 
anisotropic plasma. As discussed in Ref. this reflects the fact that such as system is 

10 



not in equilibrium and hence is unstable. The instability of a particular mode depends on 
the detailed form of the momentum distibution /o(p). It is not possible to draw general 
conclusions about the pattern of screening masses except in cases where the combination of 
derivatives 

PpT, PzTT~ 31 

op z dp p 

'his covers many smooth anisotropic distributions, in- 



has the same sign for all momenta, 
eluding the examples studied in Ref. 



1(3 ] . A negative value for this combination corresponds 
to an oblate momentum distribution with (p 2 z ) < \(p 2 p ), and a positive value to a prolate 
distribution. If the combination of derivatives has a definite sign, the integrals in H 2 = Tl yy 
can be rewritten in a form which shows that they too have the same sign. For an oblate 
distribution we get II2 < 0, indicating antiscreening while for a prolate one we get II2 > 
and screening. 

In the oblate case we also find that n o — II 3 can be written in a form which shows that 
it is negative definite. This implies that n + > and hence there is always screening for one 
component of the field. In addition, at k z = we find that n o + n 3 is negative and n 3 is 
positive, and hence n_ = II3 is positive. Since at k p = we have n_ = II3 = II 2 < 0, this 
implies that n_ must change from antiscreening to screening as the orientation of k moves 
away from the collision axis. 

For prolate distributions, in contrast, the sign of n o — n 3 is not determined and so it is 
not possible to make such definite statements about Tl±. Nonetheless the nonderivative term 
in n o (which is responsible for the usual Debye screening) is always negative and so n o 
and Iloo — S3 remain negative unless the distribution is strongly prolate. Hence we expect 
the eigenvalue n + to correspond to screening for all except the the most extreme prolate 
anisotropies. At k z = we again find that n_ = II3, but now this must be negative. At 
k p = we have n 3 = H 2 > and so both n + and H_ are positive for orientations of k close 
to the axes of distributions that are not too strongly polate. However, as the deformation 
increases, H o becomes less negative and so the range of orientations for which fl_ is positive 
shrinks. For strong enough deformations, n_ can become negative for all orientations. 

The existence of at least one unstable component of the field, for some orientations of 
k, agrees with the general results of Ref. Q]. All of the general features of II2 and Tl± 
discussed above can be seen in the numerical examples studied in Section V of Ref. 1 10], and 

n 

also in the analyses of ideally planar or linear distributions in Ref. 
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IV. CONCLUSION 



In this paper we have used a transport-theory approach to derive the polarization tensor 
of a gauge boson in a plasma which is out of thermal equilibrium. We find that the magnetic 
screening mass at lowest order is non-zero as long as the single-particle distribution function 
is anisotropic, in contrast to the more familiar case of a plasma in equilibrium. This confirms 
results previously found using a closed-time-path-integral approach. 

The full propagator for static magnetic fields in such a medium has a complicated tensor 
structure and its eigenvalues need not be positive. There can thus be antiscreening rather 
than screening for some components of the field. We have considered in detail the case 
of a plasma with a cylindrically symmetric momentum distribution. Such a distribution is 
expected to be relevant to relativistic heavy-ion collisions before thermal equilibrium has 
been reached, where the axis of collision can provide a symmetry axis for the anisotropy. In 
this case the spatial part of the static polarization tensor is real and has three orthogonal 
principal axes: k, lying along the direction the momentum carried by the gauge boson, 
y, orthogonal to the plane of k and the collision axis, and k x y. The last two of these 
have nonzero screening in general. In additional there can be an imaginary off-diagonal 
component mixing the k x y and time directions. In the special case where the momentum 
k is perpendicular to the collision axis, there is no screening for the component in the 
y-direction. 

We have also been able to derive various conditions on the signs of two of the eigenvalues 
of the tensor, under the assumption that p p ^- ~Vz^ has a definite sign. In particular we 
find that fields in the ^-direction are screened if the anisotropy is prolate, but antiscreened 
if it is oblate. We also find that there is one component which is always screened unless the 
distribution is extremely prolate. 

These differences in the screening of interactions in different directions could have impor- 
tant effects in the equilibration of the matter produced in relativistic heavy-ion collisions. 
Although only the screening or antiscreening of static fields has been examined in detail 
here, it will also be very interesting to explore dynamical aspects, such as damping rates 
and unstable modes, in these anisotropic systems 



ynamical aspect 
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